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Abstract 

When  a  laser  beam  is  directed  on  an  exciton-phonon  coupling  system 
inside  an  optical  cavity,  the  cavity  field  intensity  is  investigated  as  a 
function  of  the  driving  field  intensity  for  various  distances  between  the 
sample  and  a  metallic  surface  representing  one  of  the  mirrors  of  the  cavity. 
Optical  bistability  is  found  to  be  induced  by  the  surface  when  the  cavity 
field  is  on  resonance  with  the  excitonic  transition  frequency.  It  is  also 
found  that  at  the  threshold,  the  direction  of  switch  of  the  optical 
bistability  can  be  controlled  by  adjusting  the  distance  from  the  surface. 
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The  investigation  of  optical  properties  of  molecular  systems  near  a 
solid  surface  has  been  of  great  interest  from  both  the  fundamental  and 
practical  points  of  view.  The  study  of  an  atom  or  a  molecule  adsorbed  on  a 
solid  surface  has  attracted  much  attention^'^  in  the  past,  both  theoretically 
and  experimentally.  Exciton-phonon  coupling  systems  such  as  organic 
semiconductors  have  been  considered^  very  recently  instead  of  single  molecular 
systems  located  near  an  ideal  metal  surface.  Surf ace  -  induced  optical 
bistability  has  been  discovered  and  its  nature  and  origin  investigated. 

In  Ref.  5  we  considered,  as  an  example,  a  linear  chain  of 
polydiacetylene-toluene-sulfonate  (PTS)  oriented  parallel  to  the  metal 
surface.  The  system  closely  resembles  the  experimental  situation^  in  which  a 
single  crystal  thin  film  of  polymer  is  formed  on  the  solid  substrate. 
Polydiacetylenes  are  important  materials  for  study  in  view  of  a  variety  of 
applications  in  optical  signal  processing  such  as  ultrafast  all-optical 
switches  and  logic  gates.  One  can  foresee  potential  applications  in 
u 1 1 r ah  i  gh  -  s  p  e e d  (-TH^)  mu 1 1 ip lexing - demul t ip  1 ex ing  operation  for 
telecommunication  devices,  because  the  ultrafast  response  ( sub -p  icoseconds ) 
and  large  nonlinearity  make  these  materials  ideal  candidates. 

On  the  other  hand,  we  have  also  found  novel  optical  bistability  mediated 
by  the  exciton-phonon  interaction  in  polymers  inside  an  optical  cavity.^  It 
is  remarkable  that  the  hysteresis  loop  may  follow  the  normal  or  reversed  cycle 
depending  on  the  frequency  of  the  incident  light.  We  now  consider  in  this 
paper  a  PTS  chain  inside  an  optical  cavity.  The  chain  is  embedded  in  a  non¬ 
absorbing  material  of  dielectric  constant  and  is  oriented  parallel  to  the 
reflecting  me>,cillic  mirrors  Ciiac  constitute  the  cavity.  The  dielectric 
function  of  the  mirror  is 
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e^Cw)  “  1  -  Wp/uCaH-iPp)  ,  (1) 

vhere  is  the  plasma  frequency  and  the  damping  rate.  Our  purpose  is  to 

investigate  how  the  optical  bistability  behavior  is  influenced  by  the  presence 

of  the  metallic  surface  of  one  of  the  mirrors  if  the  polymer  is  placed  at  a 

distance  d  smaller  than  the  excitonic  wavelength.  When  an  incident  laser  beam 

is  directed  into  the  cavity,  excitons  are  excited  in  PTS  and  act  as  the 

emitting  dipoles.  The  emitted  light  is  reflected  by  the  surface,  and  the 

reflected  light  interacts  in  turn  with  the  emitting  excitons  and  changes  their 

response  to  the  incident  light.  At  the  same  time,  the  excitons  couple  with 

phonon  modes  in  PTS,  which  play  the  role  of  mediation  in  the  photon-exciton 
8 

interaction. 

To  make  the  complicated  interacting  system  manageable,  we  model  the 

excitons,  phonon  modes  and  cavity  field  as  damped  oscillators  and  assume  a 

3 

dissipative  interaction  between  the  exciton  and  surface-reflected  field.  The 

problem  is  then  treated  by  means  of  Dekker's  quantization  procedure  for 

9 

dissipative  systems.  In  this  model,  the  system  is  described  by  the  non- 
Hermitian  Hamiltonian 


H  =■  (n-ix)a^a  +  (w  - iy  )c^c  +  )  (w.-i7.)btb.  +  )  A.c^c(bt+b.) 

X  X  /  1  1  1  1  /  L  11 


- 1«  t  lu  t 


+  ig  (a^c-c^a)  +  i(a^E  e  °  -aE  e  °  )  -  /i  a^E 


1  t  * 

where  a'  (a),  stand  ror  the  creation  (annihilation) 
operators  for  the  cavity  field,  the  i-th  phonon  mode  and  the  exciton  with 
corresponding  frequencies  0,  and  respectively,  /c,  7^  and  7^^  are  the 
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damping  races  for  the  cavity  field,  exciton  and  i - th  phonon  modes, 

respectively.  A  is  the  coupling  constant  for  the  i  -  th  phonon  modes 

interacting  with  the  exciton,  and  g  is  the  exc i ton- cavi ty  field  coupling 

constant.  The  frequency  and  amplitude  of  the  coherent  driving  field  are 

denoted  by  and  ,  respectively.  The  dipole  moment  of  the  emitting  exciton 

is  represented  by  p  =  pa  with  p  as  its  matrix  element.  E  stands  for  the 

R 

reflected  field  at  the  position  of  the  dipole;  it  is  proportional  to  p  and 

hence  can  be  written  as  E_  —  E  a,  where  E  is  a  c-number.  As  wc  shall  see 

R  r  r 

later,  E^^  is  in  general  a  complex  quantity,  so  that  the  interaction  is 
dissipative . 

With  Che  dissipative  system  characterized  by  the  canonical  variables  a^ , 

ft  9 

a,  b^,  b^,  c  and  c,  we  have  according  to  Dekker  the  quantum  Liouville 

equation 

p  -  -i[a^(a.H]p]  +  i[p(H^a^l,3]  -  i  ^  [b[,[b.H]p] 

i 

+  i  ^  [p[H^b[]  .b.]  -  i(c^(c.H]p]  +  i[p[H^c^]  .c]  .  (3) 

i 

Since  we  are  mainly  interested  in  possible  optical  multistability  in  our  study 
of  the  relation  between  intensities  of  the  cavity  field  and  input  field,  we 
may  take  a  semiclassical  approach  by  neglecting  all  quantum  fluctuations.^^ 
Thus,  instead  of  the  operators  a,  b^,  c  and  c^c,  we  are  dealing  with  their 
mean  values  a  -  <a>,  -  <b^>,  ij  -  <c>  and  n  -  <c^c>.  With  these  mean  values 

defined  aj  classical  variables,  the  corresponding  equations  of  motion  can  be 
obtained  from  Eq.  (3)  as  follows: 
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a  -  - (iA^+K)a  +  g7  + 


-  -(iw.+7.)^.  -  iA.n  , 
1  1  '^1  1 


n  -  -(iA^+y)n  -  i  ^  -  ga  , 


•  *  * 
n  -  -g(arj  +a  rj)  -  27  n  , 


where  we  have  defined  the  detunings  —  Q  -  and  ^2  ”  ^  ^  ” 

7  +  7  .  The  surface  effect  on  the  exciton  is  reflected  in  the  surface- 

X  s 

induced  frequency  shift 


Wg  -  |m|  Re(E^/fi) 


and  the  decay  rate 


7g  -  |m| 


The  reflected  electric  field  in  the  dipole  direction  can  be  found  by 
a  classical  approach.  We  take  the  oscillating  dipole  moment  to  be  located  at 
a  distance  d  froir  the  surface  of  a  semi- inf inite  metal.  From  electromagnetic 
field  theory,  we  find  in  a  straightforward  manner  that 


3  1^  2  2 

-  -iJT^  J  du  u  R|j  exp(2i^^x)/yi-u 
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In  the  derivation  of  Eq.  (6),  we  have  assumed  that  the  dipole  is  perpendicular 

to  the  surface.  Here  we  have  also  defined  k  -  w  /c ,  x  -  /7,  k  d  and 

ox  1  o 


^^2  -  ^2^1 
^  ^2^1 


^2  “  ^  ,  i  -  1,2 


Substituting  Eq.  (6)  in  Eqs.  (5),  we  find 


3  1  3 

w  -  r  Im  du  u  R||  exp(2i^,x)//i, 
s  2  J  Q  II  1  J 


^  J'®  « 

3  3 

-Re  du  u  R|  exp(2i/ij_x)/Mj^  , 


2  /—  2  3  3 

which  have  been  expressed  in  the  unit  of  7^  -  —  JT^  [mI  ■  Ftora  now  on, 
and  throughout  this  paper,  all  energies  will  be  measured  by  this  unit.  If  we 
introduce  a  dimensionless  input  field  intensity  I^^  —  g  1e^|  find 
from  Eqs.  (4),  after  some  algebraic  manipulation,  that  the  mean  exciton  number 
n  satisfies  the  cubic  equation 


Ap[l+(A^A)^]n^  +  2A  [A^(g^-A^A2)A^'A2]n^ 


+  {(1+7^  +  (g^-A  A,)A]^  +  (A  +  (1+7  )A  /«:]^)n  -  I.  -  0 
s  I  L  I  Si  in 


2  2 

At  the  same  time,  the  cavity  field  intensity  I  -  |q|  /y  can  be  expressed 

CSV  X 

in  terms  of  the  input  field  intensity  I.  as 
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(1+7)^  +  (A  -A  n)^ 

^cav  “  - 2 - 2 -  - 2  ^in  ’ 

[1+7^+g  A-A.  (A  -A  n)//c)  +  [A  (1+7  )A+Ao--^„n] 


■»/%  ‘ 
2  P 


where  we  have  defined  A 


-  2  I 


In  the  limit  /c  -►  ®,  we  are  dealing  with 


a  bad  cavity,  namely,  the  cavity  field  damps  out  almost  immediately. 
Equations  (9)  and  (10)  then  reduce  to 


A^n^  -  2A  A_n^  +  ((I+7  +  A^jn  -  I. 

p  p  2  '3  2 '  in 


(11) 


cav  — 


I. 

in 


(12) 


which  are  identical  to  the  results  of  Ref.  5,  as  expected. 

Equation  (9)  is  solved  numerically,  and  the  solutions  are  used  to  study 

the  tristability  of  the  cavity  field  intensity  in  Eq .  (10)  under  various 

conditions.  When  the  incident  field  frequency  is  on  resonance  with  both  the 

cavity  field  and  the  excitonic  transition,  the  occurence  of  optical 

bistability  d.^pends  on  the  location  of  the  chain  in  the  cavity.  If  the  chain 

is  far  away  from  the  mirror  surfaces  of  the  cavity,  Eq .  (9)  has  only  one 

physical  solution  and  there  is  no  optical  bistability.  If  the  PTS  chain  is 

near  a  mirror  surface  or  if  d  is  smaller  than  the  excitonic  wavelength,  Eq . 

(9)  may  possess  three  real  solutions,  and  it  is  possible  to  find  optical 

bistability  in  the  cavity.  As  discussed  in  Ref.  5,  this  surface- induced 

bistability  is  caused  by  both  the  frequency  shift  w  and  the  width  change  7 

s  °  s 

of  the  excitonic  transition  due  to  the  surface.  If  the  distance  becomes  too 
small  or  X  «  1,  the  surface  absorption  becomes  so  strong  that  the  lifetime  of 
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Che  exciton  is  too  shore  to  interact  strongly  and  nonlinearly  with  the  phonon 
modes.  Consequently,  there  can  no  longer  be  bistability. 

Off  resonance,  the  behavior  depends  very  much  on  the  sign  of  the 

detuning.  When  >  0,  Che  optical  bistability  is  expressed  in  a  normal 
5 

hysteresis  loop,  and  the  presence  of  a  surface  can  not  change  this  nature. 
However,  a  nearby  surface  can  affect  the  threshold  input  intensity  and  the 
contrast  of  the  hysteresis  loop.  When  <  0,  optical  bistability  can  usually 
be  found  with  the  inverted  hysteresis  loop,  even  if  the  sample  is  far  away 
from  the  surface.  In  the  presence  of  a  nearby  surface,  our  calculations  show 
a  completely  new  phenomenon,  namely  that  the  two  cavity  field  intensities 
cross  each  other  as  the  driving  field  intensity  changes. 

In  our  numerical  work,  we  have  included  all  four  phonon  -  .modes  that 

7  8  11 

couple  most  strongly  to  the  excitons.  ’  ’  These  i.’.odes  are  characterized  by 
the  parameters  -  5.16,  -  2;  -  3.68,  -  1.66;  -  2.98,  A^  -  0.46; 

--  2.36,  A^  -  0.48.  Some  of  the  results  for  A,  <  0  are  plotted  in  Fig.  1  in 
which  we  demonstrate  the  cht^nge  of  th--  bistability  nature  by  varying  the 
distance  from  the  surface  only.  Figure  1(a)  shows  no  bistability  because  of 
the  small  distance.  As  mentioned  above,  the  strong  absorption  by  the  surface 
prevents  the  exciton  from  interacting  with  phonon  modes  nonlinearly  There 
can  not  be  opt?'’al  bistability  for  any  A^^.  In  Figs.  1(b)  and  1(c),  we  find 
the  interesting  new  phenomenon  where  the  two  output  states  cross  each  other, 
forming  a  peculiar  8-shaped  loop.  The  cavity  field  changes  the  normal 
bistability  as  the  driving  field  intensity  increases,  but  '^ollows  the  inverted 
loop  as  Ij^^  decreases.  The  hysteresis  loop  becomes  a  triangle  in  Fig.  1(d)  at 
the  distance  x  -  2.0.  After  this  point,  the  loop  becomes  inverted  and  stays 
that  way  as  the  distance  increases  further. 
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As  discussed  in  Ref.  7  in  more  detail,  in  an  optica^  cavity,  the  real 
part  of  the  optical  response  of  PTS  or  the  dispersion  part  plays  a  more 
important  role  than  the  absorption  part  because  the  refractive  index  is 
changed  by  the  excittns  which  are  excited  by  the  j-npu*"  ^ield.  Our  numerical 
computations  indicate  that  the  major  contribution  to  the  bistability  comes 
from  the  exciton  energy  shift  induced  by  the  surface.  This  is  in  contrast  to 
the  case  discussed  in  Ref.  5,  where  the  absorptive  bistability  is  caused  by 
reduced  vacuum  fluctuations. 

This  research  was  supported  in  part  by  the  Office  of  Naval  Research  and 
the  National  Science  Foundation  under  Grant  No.  CHE-9016789. 
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Figure  Captions 

1.  Cavity  field  intensity  I  versus  the  incident  field  intensity  I.  for 

cav  in 

various  distances  x  between  this  system  and  the  metal  surface  when  “ 
-0.25,  u  /w  -  0.8,  A.  -  6,  g^/k  -  0.5.  (a)  x  -  0.5,  (b)  x  -  0.75,  (c)  x 

X  p  2 

-  1.25,  (d)  X  -  2.0,  (e)  X  -  2.5,  (f)  x  -  ®. 


